Measuring entanglement using quantum quenches 
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We show that block entanglement entropies in one-dimensional systems close to a quantum critical 
point can in principle be measured in terms of the population of low-lying energy levels following a 
certain type of local quantum quench. 



Entanglement, and its quantification through entangle- 
ment entropy, have become increasingly important tools 
in the study of quantum many-body systems. The fact 
that the entanglement entropy of a region of linear size 
£ in the ground state of a system with short-range inter- 
actions in d dimensions grows like £'^~^ [1] (up to possi- 
ble logarithms), as compared with a typical state where 
it is extensive, both explains the success of modern nu- 
merical methods and informs their further development. 
Moreover it gives a basis-independent way of detecting 
quantum critical behavior and topological phases. 

Briefly, the entanglement entropy is defined as fol- 
lows: given a bipartition of the Hilbert space Ti = 
Ha ®T~Lb (which usually corresponds to the degrees of 
freedom lying in mutually exclusive regions A and B of 
rf-dimcnsional space), the reduced density matrix of (say) 
A is given by pA — Tr^^^ |0)(0|. The Renyi entropies are 
then given by 5^"' = (1 — n)^Mog Tr-H_4 p^. The von 
Neumann entropy — Tr^^ pA log pA is formally the limit 
as n — > 1. For different n these encode the entanglement 
spectrum of pA ■ if we make a Schmidt decomposition 



|0) = 



(1) 



then the non-zero eigenvalues of p^ are given by cf.. On 
the other hand Tr-u^ Pa — J2k "^fe"- state has low 

entanglement then only a few eigenvalues are apprecia- 
ble and the entropies are small, while if the state is max- 
imally entangled the entropies are of the order of the 
lesser of the dimensions of Ha and H b ■ 

In exactly solvable models, and more generally for sys- 
tems at or near a quantum critical point, considerable 
progress has been made in relating the behavior of the 
entanglement entropy to other universal data of the un- 
derlying long-distance theory, typically a quantum field 
theory (QFT) 0. In particular, for d = 1, the coefficient 
of the 'area' law (in this case the number of boundary 
points between A and B) is predicted to diverge as log^, 
where ^ is the correlation length, as the critical point is 
approached. At the critical point, the entropy of a finite 
interval of length £ immersed in a much larger system 
grows like log£. When the low energy, long wavelength 
physics is described by a conformal field theory (CFT), 



the coefficient is given by the central charge, or confor- 
mal anomaly number, c: the Renyi entropies behave as 
d 13 S'x^ - (c/6)(l + n-^)log£. Many other predic- 
tions have been made by now. For example, when A 
consists of two disjoint intervals, the entropy encodes all 
the data of the CFT, which is therefore recoverable solely 
from properties of the ground state wave function [sl, Q . 
These predictions, both in d — 1 and higher, have been 
well verified in numerical investigations, and indeed the 
low degree of entanglement is partly responsible for the 
success of DMRG and tensor network methods 't'I . 

However it seems difficult to devise a method whereby 
entanglement entropy in an extended system of the type 
discussed above could be measured, even in principle, in a 
real experiment, since it is intrinsically non-local. While 
various methods have been proposed for measuring en- 
tanglement in systems with a finite dimensional Hilbert 
space 0, their complexity increases with the system size. 
What we have in mind is a conventional condensed mat- 
ter experiment whose difficulty does not in principle in- 
crease with the system size. In some simple systems the 
entanglement entropy can be indirectly recovered from 
an (in principle) measurable correlation function [qI. Iiol 
or the distribution of suitably chosen observables |ll| . 
but this connection is system-specific. 

Various suggestions have been made as to how charge 
H [11] or number fu] lluctuations in the subsystem 
A, either in space or in time, may provide a measure- 
ment of entanglement entropy. However, the generality 
of these observations has been questioned [31 ■ They are 
restricted to systems with a conserved current, and it 
is difficult to see how these ideas could measure entan- 
glement of neutral degrees of freedom, or apply in cases 
when there is no such conservation law. More seriously, 
the logarithmic behavior which appears in these analyses 
can be traced to the fact that, in one dimension, current- 
current correlators {J(xi)J{x2)) behave like \xi — X2\~'^ 
on separations much smaller than the size of the subsys- 
tem (but larger than the microscopic scale, and similarly 
in the time domain), giving rise to logarithms on inte- 
gration over xi and X2. Although the functional form of 
such logarithms turns out to be identical to that of the 
entanglement entropy for simple geometries, the analysis 



of Ref. ^ shows that, even in the slightly more compli- 
cated case where A consists of two disjoint intervals, the 
entanglement entropy has a form which cannot be sim- 
ply expressed as an integral over correlation functions of 
local fields. 

The proposal in Ref. [12 is an example of a local quan- 
tum quench T^, l3| , whereby the hamiltonian of a quan- 
tum system is instantaneously changed H ^ H' in a. local 
way, so that the quantum state, which was the ground 
state of H, now evolves according to H' . In general, the 
real time behavior after such a local quench is relatively 
simple [17| : the additional energy near the quench site is 
radiated away as the quasiparticles of H' , moving semi- 
classically at their group velocities Vg . These also propa- 
gate changes in entanglement through the system, which, 
in a critical Id system, therefore grow like \og{vgt). 

In this letter we propose a different type of local quan- 
tum quench which, in principle, directly measures the 
Renyi entropies. We show that these are simply given by 
the probability Pq = Iff' (0|0)hP of finding the system in 
the ground state of H' after the quench. This quantity 
is called the fidelity, and it has been used extensively in 
characterizing quantum critical behavior (see, for exam- 
ple, Ref. [lih . but it has largely been restricted to global, 
rather than local changes in the hamiltonian. (See, how- 
ever, [l9j.) It may be objected that Pq is not directly ex- 
perimentally observable, even in a gapped system. How- 
ever we also show that, in systems close to a quantum 
critical point described by a QFT with a linear disper- 
sion relation, the probability P{E) of finding the system 
in a low-lying excited state of energy E is given by Pq 
times a calculable factor. In principle, P{E) is measur- 
able if the system is coupled weakly in a known fashion 
to other modes whose spectrum can be analyzed. This 
then gives an indirect measurement of Pq and hence the 
Renyi entropy. The analogof P{E) for global quenches 
has been discussed in 2^ 2ll . 



Let us define more precisely the local quench. In com- 
puting Renyi entropies, it is convenient to imagine n 
copies of the original system. We now suppose that these 
copies can actually be manufactured, to a tolerance to be 
discussed later. Consider therefore n identical copies of 
the system so that the fuU Hilbert space is "H = (8>"^i'Hj. 
Initially they are decoupled so that the hamiltonian is 
H = ^^'^ each system is in its ground state (as- 

sumed unique), so the combined system is in its ground 
state \Q)h — Yij Consider the same bipartite de- 

composition of each Hj = HjA <8) 'Hjb- All the Hja are 
isomorphic. Let n„ be the permutation operator (uni- 
tary on H) which maps Hja "Hq+i)^ (mod n) and 
acts as the identity on all the T-ljB- (Such an operator 
for the case n = 2 was introduced in Ref. [H where it is 
called a swap operator.) Then it is well-known (see e.g. 
Ref. [2^ for the case n = 2) that using the decomposition 
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FIG. 1: The action of a single twist operator on two one- 
dimensional systems with short-range interactions. 



FIG. 2: A single twist is equivalent to the closing of a point 
contact when n is even and the ground state is invariant under 
reflection. 



([T]) and the orthonormality of the Schmidt states. 



H(o|n„|o)H 
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Tr«^ Pa 



(2) 



On the other hand, if we define H' — n^^i/n„, H' and 
H are isospectral, and \Q)h' = ^n^\^)H is the ground 
state of H' . Hence, thinking oi H ^ H' as a. quantum 
quench, the modulus squared of ^ is nothing but Pq — 

\H'{mH?- 

The point about this elementary observation is that, 
when H has only short-ranged interactions, and A and 
B are spatially disjoint regions, the difference between 
H and H' is restricted to the boundary between A and 
B. Thus while the action of H„ on states, expressed in 
particular basis, may be complicated (see, e.g. Fig. 1 of 
Ref. 1221 ) , its action on the combined hamiltonian is sim- 
ple. As a example, consider two copies of a Heisenberg 
spin chain with nearest neighbor interactions and hamil- 
tonian J'Y^i (J I ■ (Ti+i. Take A to be the set of sites with 
I < and B those with / > 1. Then 



Tjl ZJ T (1) (2) I T l^J IJ-J T IJ-J IJ-J T 

H —H — Jctq -(Tl '+J<Tq -al — JCTq -al — JCTq -ct} 

This is illustrated in Fig. [TJ (For an odd number of 
such twists it is also necessary to twist the boundary 
conditions at the same time. This has no effect in open 
systems or for periodic boundary conditions in systems 
much larger than the correlation length, but is impor- 
tant for finite ungapped systems with periodic boundary 
conditions.) Such twist operators have been studied for 
integrable spin chains in i24| . 

Note that for an infinite system or one with periodic 
boundary conditions, one of the chains can be parity re- 
versed as in Fig. [2] (assuming the ground state |0) is in- 
variant under this reflection) to obtain a local quench 
corresponding to the closing of a point contact coupling 
to two external leads on either side. In a real experiment. 
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of course, it may be difficult to ensure that the new cou- 
plings added in H' are precisely the same as those deleted 
in H . However, at least in a gapped system, the results 
should be robust to such imperfections as long as they 
are on scales smaller than the gap, or the energy E of 
excited states. 

The operators n„ corresponding to a local modifica- 
tion of H at some point I in a one dimensional sys - 
tem are called twist operators Tn{l) in the literature [25|. 
From the above we see that the Renyi entanglement en- 
tropy between the left and right halves of such a sys- 
tem is given in terms of the ground state expectation 
value _ff (0|7^i(0)|0)_ff- Simfiarly the Renyi entanglement 
entropy between an interval of length i = \li — hi and 
the rest of the system is given by /f(0|7^(Zi)7^(/2)|0)//, 
and so on. In all cases, these matrix elements may be 
viewed equivalently as the fidelity following a quantum 
quench of H to H' ^ Hfe TnihV H Hfc Tnih), or equiva- 
lently the probability Pq that the system is found in the 
ground state of H' when its energy is measured. 

In one-dimensional systems close to a critical point, 
considerable effort has gone into analyzing the behavior 
of correlators of these twist operators, especially for sys- 
tems with dynamical scaling exponent z = I, which are 
described in the scaling limit by a relativistic QFT, and 
at the critical point, by a CFT m. In these cases, it 
follows from the early work of Ref. y that Tn has scaling 
dimension a;„ = (c/12)(n — where c is the conformal 
anomaly number of the CFT. By scaling, this means for 
a single twist in a system with a finite gap A, Pg scales as 
&^A^^" as A ^ 0. Similarly, for an interval of length £ at 
the critical points, it behaves like &„€~^^". The constants 
bn and b'^ are not universal, but their ratio is expected 
to be, in units where v = 1 [28j. 

However, measuring the total energy of a many-body 
system to accuracy 0(1) is unfeasible even when the 
spectrum is gapped. In fact, after such a local quench, 
most of the excess energy goes into states whose energy 
is of the order of the inverse cut-off, or band width. 
This can be seen by studying the time dependence of 
the energy density, given by the component Too(a;,t) of 
the energy momentum tensor, following the action of a 
single twist. This is given in the Schrodinger picture 
by //(0|r„(0)te*^*Too(x)e-*^*r„(0)|0)/f. Onjhe other 
hand, in a CFT, we can write Too = T + T, and, in 
imaginary time 

ff(o|r„(o)te-^-'r(x)e-^-r„(o)|o)^f 

^ (e/24)(l-n-2)(T + T0^ 
{x + iT)^(x — ir')^ 

This follows from Eq. (11) of 4J. A similar result holds for 
(T). Continuing this result to r = it + S and r' = —it + 6, 
where (5 is a UV cut-off of the order of the lattice spacing 



or inverse band width, we see that, after the quench, 
(Too(x,t)) - -^(n-n-'^) ( ^ + {x -x}\ 

in units where v = 1. The total energy j {Too)dx diverges 
like S~^, and, in a CFT, it is concentrated in a region of 
width 0{S) near the light cone. In a more general lattice 
theory we expect the energy to be carried off by quasi- 
particles moving semi-classically with their appropriate 
group velocities. Similar considerations apply to changes 
in the entanglement entropy 

We now argue that, despite the fact that most of the 
energy is radiated in this non-universal manner, the pop- 
ulation P{E) of the states with energies E much less 
than the cut-off is universal, and is directly related to 
the Renyi entropy. This because its Laplace transform is 
given by correlation functions of twist operators in imag- 
inary time: 

J P{E)e-^-dE = H{0\e-^"'-"^lO)H 

= h{0\{1[ Tnih, tY) ( n Tn^fc, 0)|0)h 
k k 

We now consider various simple cases of this. 

Single twist in an infinite ungapped system. At the 
critical point the two-point function (7^(0, r)'''7^(0, 0)) ^ 
6„r~^^", where is the same constant appearing in 
the result for the Renyi entropy of an interval of length 
£ (in units where v — 1.) In this case the spec- 
trum of H' is continuous and we see that P{E)dE ~ 
{bn/T{2xn))E'^'="~^dE. Note that the probability of 
finding the system in a state of energy < E approaches 
zero as £■ — >■ 0. This is an example of an 'orthogonal- 
ity catastrophe' as was first observed in the X-ray edge 
singularity [26| . As in that problem, the probability of 
occupation of low-lying states obeys a universal power 
law [27| . 

Single twist in a gapped system. When the theory is 
gapped, the 2-point function behaves asymptotically like 
(Tn)]"^, giving a probability Pq ~ 6^ A^^" that the system 
is in the ground state. More interesting is the occupation 
of low lying states. In [gl] it was shown that Tn couples 
to multi-particle states in the continuum. As an example, 
the coupling to 2-particle states with E > 2A gives 

n n 

P{E) ^PoJ2 EE l^^^(^i'^2)r<5(i? - E,, - Eu,) 

ki ,^2 i— 1 i — 1 

where Py is a form factor coupling 7^ to a particle of 
momentum ki in copy i and of momentum k2 in copy 
j. For integrable models these are calculable Q and 
thus measurement of P(P) gives direct access to Pq and 
therefore the Renyi entropy. 

Single twist in a finite ungapped system. Suppose 
the initial system has finite length L with open bound- 
ary conditions. The form of the 2-point function 



(7^(a;, r)'''7^(a;, 0)) can be found by a conformal map- 
ping to the upper half plane. As r — > oo we find 
^0 — &n.9ri ((^/tt) sin(7ra;/L))~^^", where x is the dis- 
tance of the twist from one end of the system, and log 
gives the boundary entropy. This agrees with the result 
in Q for the Renyi entropy. The coupling to the excited 
states can be found using the methods described in Q. 
The leading corrections as r — >■ oo come from two exci- 
tations each of energy irXg/L (where boundary 
scaling dimension) propagating in two different copies of 
the available n. We then find that 

P{E = 2t:Xs/L) = Podnix,) ((L/^) sin(^x/L))^"'' 

where dn{xs) is a known function of j;s Q. 

Multiple twists. This is more complicated still, since 
for p twists we need to know the 2p-point function of 
twist fields. However analytic results are available from 
CFT for various limiting cases. For example take the 
case p = 2, that is an interval A of length £ in an infinite 
system. As discussed above, the probability Pq of finding 
the system in its ground state gives the Renyi entropy. 
For i <^ T, that is energies E ^ the methods of 
show that the product of two twist operators can be 
written as an infinite sum over all possible local scaling 
fields of the CFT on each of the n copies. The leading 
correction comes from when two of these have dimension 
X and all the rest correspond to the identity. Thus 
we have 

P{E)dE = PqY, dn{x){tEf''-^dE + 



where c?„ [x) is once again known, the sum is over all the 
bulk scaling dimensions, and the neglected terms, of or- 
der (^i<;)2(a;i+a;3+a:t) g^j.g ^^jg^ calculable and, for n > 2, 

encode further CFT data. (For n = 2 a closed form for 
P{E) is available since the 4-pt function of twist oper- 
ators is related to the partition function of the CFT on 
a torus 0.) For E ^ we find the result for single, 
independent twists. 

Another important example of a local quench, however 
not related to the Renyi entropy, is where two pieces of a 
system are joined together to form a larger one. The real 
time evolution after such a quench was studied in (T"" 



in 



More recently the fidelity was evaluated using CFT meth- 
ods [l^l ■ We can also use evolution in imaginary time to 
predict the population of low energy states following such 
a quench. In this case, instead of a twist operator, we 
have a slit in the world sheet, whose scaling dimension 
in a CFT was computed in [29| to be c/16. Thus, in an 
infinite system, we have P{E)dE cx E'^'^/^'^-^dE. 

Although the purpose of this article has been to show 
only that it is possible in principle to measure Renyi en- 
tropies, it is interesting to consider whether this proposal 
is at all practical. As mentioned earlier, imperfections in 
the faithfulness of the copies or of the precise details of 



the quench should be unimportant as long as they occur 
at energy scales much less than the gap A or the energy E 
above the ground state. In principle, P[E) can be mea- 
sured if the system couples weakly, in a known manner, 
to other degrees of freedom whose energy distribution 
can be spectrally analyzed. For example, if the excited 
states of H' can decay through a coupling to the elec- 
tromagnetic field, P{E) should be recoverable from the 
photon spectrum. In optical lattices, the energy E could 
be transferred to kinetic energy of the atoms, whose spec- 
trum could then be analyzed after removal of the trap. In 
both cases, in order to produce an observable spectrum, 
it would of course be necessary to consider a small but 
finite density of twists, produced continuously at a low 
rate. Under such conditions, the above analysis shows 
that we should obtain simple universal behavior, charac- 
teristic of independent twists, for energies E larger than 
mean density (in units where v — 1), but much smaller 
than the band width. 

Although we have restricted attention to case of one 
dimension, where twists are local, the analysis in prin- 
ciple extends to higher dimensions. In particular it can 
be seen that the action of twist operators can change 
the topology and so reveal the entanglement entropy of 
topological phases. We have throughout assumed that 
the ground state |0) is unique, and there are interesting 
consequences, even in one dimension, when this is not 
the case. 
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